Scattering of Solitary Waves in Granular Media 
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A detailed numerical study of the scattering of solitary waves by a barrier, in a granular media with 
Hertzian contact, shows the existence of secondary multipulse structures generated at the interface 
of two " sonic vacua" , which have a similar structure as the one previously found by Nesterenko and 
coworkers. 

PACS numbers: 46.40.Cd; 45.70.-n; 47.20. Ky 



INTRODUCTION 

Nesterenko noticed that the propagation of a per- 
turbation in a chain of beads in Hertzian contact pos- 
sesses soliton-like features. Several studies, theoreti- 
cal as well as experimental, 0, Q 0], have con- 
firmed the existence of such soliton-like pulses. De- 
spite the great deal of recent work on the subject 
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granular media remains a challenge and new effects are 
there to be discovered and studied. 

The simplest granular systems are one-dimensional 
chain of elastic spheres. It is known that in the 
regime where Hertz's static solution for the contact of 
elastic spheres applies, the spheres may be considered 
as point masses interacting through massless nonlinear 
springs with elastic force F = kS 3 ^ 2 , where S is the over- 
lap of contacts and k is the spring constant (a function 
of the material properties). Let Vi{t) represents the dis- 
placement of the center of the i-th sphere from its initial 
equilibrium position, and assume that the i-th sphere, 
of mass m, has neighbor of different radius (and/or me- 
chanical properties). Then, in absence of load and in a 
frictionless medium, the equation of motion for the i-th 
sphere reads 



FIG. 1: 

NUMERICAL STUDY AND RESULTS 

Consider a set of spheres with two different radius o 
and b. It is known that adjacent spheres of radius a and 
6 will interact with a force F — kS 3 ^ 2 , where 
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where it is understood that the brackets take the argu- 
ment value if they are positive and zero otherwise, ensur- 
ing that the spheres interact only when in contact. 

Nesterenko, Lazaridi and Sibiryakov [l6| studied ex- 
perimentally as well numerically the interaction of a soli- 
tary wave with the boundary of two " sonic vacua" using 
a setup similar as the one shown in Figure 1, with a wall 
at the right hand side. They also analyzed the case where 
the small and big beads are interchanged. 

In this work we make a detailed numerical study of the 
propagation of solitary waves in a linear chain of beads 
composed of two types of bead sizes, as shown in Figure 
1, repeating the analysis done in 16]. It will be assumed 
that all spheres are made of the same material and that 
both ends of the chain are free to move. 



and E is the Young modulus and v the Poisson ratio of 
the bead material. 

We will consider the scattering of solitary waves in a 
setup as that shown in Fig. 1, similar to the one used in 
^f|[2(]], where there are M beads, of which N has radius 
a and mass mi and L radius b and mass m 2 ; a < b. 
The displacements of the beads are governed by a set 
of eqs. of motion that can be readily obtained from the 
successive application of eq. JTJ , having in mind that the 
cq. of motion for the first (resp. the last) sphere only 
includes the second (resp. the first) term, in case there 
is no wall (as we here assume). 

For numerical convenience, we rescale the variables, 
defining 



mi 



(4) 



In the first part of our analysis of the effect of the in- 
terface on the dynamics of solitary waves, the parameters 
of the system are chosen as N = 60, L = 60, a = 1, b = 2, 
TOi/m.2 = 1/8. We assume that initially all beads are at 
rest, except for the first bead at the left side of the chain. 
This bead is supposed to have a nonzero value of velocity 
in order to generate the soliton-like perturbation in the 
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chain. That is, we assume that the initial conditions are 

Ui (0) = 0, i=l,...,M, ui(0) = 1 
u x (0) = 0, i = 2,...,M, 

In this case, the scattering process is as follows: the 
solitary wave arrives at the interface between the two 
"sonic vacua" and it splits into a transmitted and a re- 
flected solitary wave as already observed in . In addi- 
tion to what reported there, we have observed the pres- 
ence of extra multipulse structures. 
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FIG. 2: Multipulse structure produced by the scattered soli- 
tary wave (by a setup as that shown in Fig. 1) as a function of 
the bead number, at different times. The vertical line denotes 
the interface position. 

In Figure 2 we have plotted the velocity of the pulses 
as a function of the bead number for different times 
t = 140, 165, 190, 220 and plots where shifted by 0.025 
for clarity. At t — 140 it is observed (part of) the re- 
flected solitary wave and behind it two pulses appear. In 
addition, an unresolved perturbation around the inter- 
face is present. As shown in the figure, at later times 
multipulse structures appears. 

The perturbations in the multipulse structure possess 
much less energy than the solitary waves; comparing the 
velocity of the reflected solitary wave with the velocity 
of the first pulse in the multipulse structure shows that 
in the case at hand, the perturbation has around 3% of 
the velocity of the solitary wave. 
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FIG. 3: Velocity of bead 54 as a function of time. 



Figure 3 shows the velocity of bead 54 as a function of 
time. The large peak corresponds to the reflected solitary 
wave, while the velocity of the incident solitary wave is 
not shown. The two multipulse structures observed here 
have they counterpart in the structures observed at t — 
240 in Figure 2. 

It is also interesting to study the dependence of the ef- 
fects that we have observed on the strength of the impact 
by the first bead. To that end we use as initial conditions 

Ui (0) = 0, i=l,...,M, u x (0) = 3 
u x (0) = 0, i = 2,...,M. 
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FIG. 4: Velocity of bead 54 as a function of time. 

Figure 4 shows the dependence of velocity of beads 54 
as a function of time, for these initial conditions. As 
expected, all perturbations are enhanced and because of 
the larger magnitude of the initial velocity, perturbations 
are shifted to earlier times as compared with Figure 3. No 
big differences are found between both cases, indicating 
that (at least to these orders of magnitude) the impact 
velocity has no incidence in the the appearance of the 
multipulse structures. 

Before continuing with our a analysis we remark that 
since the effects that we observe are small it is impor- 
tant to be sure of the accuracy of the numerics. We have 
numerically studied the system of equations by using an 
explicit Runge-Kutta method of 5th order based on the 
Dormand-Prince coefficients, with local extrapolation. 
As step size controller we have used the proportional- 
integral step control, which gives a smooth step size se- 
quence. We have tested our numerics by comparing with 
results given in the literature as the one shown in |6J and 
[l7j . obtaining a good agreement with them. In addi- 
tion, as a crude check of numerical integration accuracy 
we have found that the final kinetic energy of the trans- 
mitted solitary wave in the heavy system differs from its 
initial kinetic energy in that system by 10~ 7 . 

When studying the dependence of the multipulse struc- 
ture as a function of the ratio between the beads, one ob- 
serves that the number of pulses in the multipulse struc- 
ture increases when the mass ratio m x /m2 decreases and 
viceversa. Also notice that for lower mass ratios m\jmi 
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the multipulse structure take more time to appear. This 
is shown in Figures 5 where the velocity of bead 54 as 
a function of time is plotted. Mass ratios are chosen as 
™i/ TO 2 = 1/1-6 3 , m 1 /ni2 = 1/2 3 and mi/m 2 = 1/2. 5 3 . 
The curves corresponding to the last two mass ratios have 
been shifted by 0.05 units for clarity. 
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FIG. 5: Velocity of bead 54 as a function of time for three 
different mass ratio. 

One could ask whether similar structures also appear in 
the case when the solitary wave moves from the " heavy" 
system to the "light" system (from right to left in Figure 
1). This system was first analyzed in and here we 
reanalyze their numerical calculation. To this end we 
choose the system's parameters as a = 2, b = 1, mi/m-2 = 
8, with 30 big beads and 170 small ones and the initial 
conditions J5J|. 



it appears a second multipulse structure, with similar 
characteristics than the first one but with less energy. 
This is shown in Figure 7. 
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FIG. 7: Grain displacements with respect to their original 
equilibrium position and pulse velocity as a function of bead 
number at t — 115. The pulse velocity has been scaled by a 
factor 6. 

Notice that there are some similarities between this 
effect and the one observed before. Indeed, when the 
perturbation travels from the light to the heavy system, 
the multipulse structures are generated at the interface, 
leaving some energy behind the interface. If we zoom-in 
Figure 7 we notice that something similar happens here. 
Indeed, in Figure 8 and 9 we show how the secondary 
multipulse structure is generated by a "seed" of energy 
that remains at the interface after the primary multipulse 
structure has been generated. 
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FIG. 6: Snapshots at times ranging from t = 34.4 to t — 40 
showing the way the pulse velocity deforms to produce the 
first multipulse structure. Plots were shifted by 10 units for 
display reasons. 

The scattering process is as follows: the solitary wave 
comes to the interface and no backscattered solitary wave 
is observed but a multipulse structure is generated as 
shown in Figure 6. This is a remarkable phenomenon 
whose origin is not completely clear to us (apart from 
the fact that it originates in the discreteness and nonlin- 
earity of the system in question). It was first observed 
by Nesterenko and coworkers in [16j. 

The new effect that we report here is that a time after 
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FIG. 8: Zoom-in of Fig. 7, at time t = 60, showing the 
"seed" of energy behind the behind the interface and part of 
the primary multipulse structure. 

Using a detailed numerical approach, we have stud- 
ied the scattering of solitary waves in a system with to 
"sonic vacua" in a granular media with Hertzian contact. 
We have studied the same system as the one analyzed in 
(ljj (except for the existence of a wall, that makes no 
difference in the final results) and found that, in addi- 
tion to what observed there, multipulse structures with 
smaller energies emerge. As far as we know, this kind of 
secondary multipulse structures have not been observed, 
yet, in experiments. We are aware of a new experimental 
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FIG. 9: Zoom-in of Fig. 7, at time t = 88, showing the 
generation of the secondary multipulse structure. 



method suitable for the investigation of solitary waves at 
walls (we thank Prof. Francisco Melo for giving us a copy 
of his paper after the completion of this work) . We 
expect that this method could be used for studying the 
effects shown here. 

Although we have given some information about the 
formation of the secondary multipulses structures shown 
here, their origin is not completely clear to us but it is 
interesting to notice that similar multipulse structures 
appear in the scattering of two solitary waves d (al- 
though in that case only primary multipulse structures 
appear). We agree with them that part of the explana- 
tion lies in the discreteness of the medium through which 
the solitary waves propagate. 

Finally, it is worth to mention that although the effects 
shown here are small, in our opinion they are important 
in the sense that they form part of the dynamics of the 
scattering of solitary waves by interfaces and then deserve 
to be studied and their origin further clarified. 
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